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CARRYING THE MASS FLUX TERMS EXACTLY 
IN THE FIRST AND SECOND MOMENT EQUATIONS 
OF COMPRESSIBLE TURBULENCE 


J.R. Ristorcelli 1 

Institute for Computer Applications in Science and Engineering 
NASA Langley Research Center 
Hampton, VA 23681 


ABSTRACT 

In compressible turbulence models it is assumed that the Favre-mean velocities are suit- 
able approximations to the Reynolds-mean velocities in order to close unknown terms. This 
neglects, in the mean momentum and energy equations, the contribution to the stress and 
work terms by the mean of the fluctuating Favre velocity, a quantity proportional to the 
turbulent mass flux. As the stress and work terms do not introduce any new unknown cor- 
relations requiring closure in either k — e or Reynolds stress closures and because the exact 
form of the terms can, with little additional work, be carried there is no need to make any 
modeling assumptions. In the Reynolds stress equations the viscous terms appear naturally 
in Reynolds variables while the problem is posed in Favre variables. In the process of split- 
ting the viscous terms into the viscous transport terms, carried in Favre variables, and the 
dissipation terms, carried in Reynolds variables, important contributions from the mass flux 
appear. The accurate accounting of these terms is important for any consistent near wall 
modeling and the retention of the mass flux terms is important in complex compressible 
turbulent flows. 


1 This research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NAS1-19480 while the author was in residence at the Institute for Computer Applications in 
Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23681. 
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Introduction 

The stress and work terms in the mean momentum, mean energy and the Reynolds 
stress equations are usually modeled by assuming that the Favre mean velocity is a suitable 
approximation to the Reynolds mean velocity. This neglects the contribution of the mean 
of the fluctuating Favre velocity which, related to the turbulent mass flux, quantifies the 
difference between the Reynolds and Favre mean velocities of the mean flow. As the stress 
and work terms do not introduce any new unknown quantities in k — e or Reynolds stress 
models (an eddy viscosity expression for the mass flux is typically carried) and because 
they can be carried exactly with little additional complexity there is no need or justification 
for any modeling assumptions. The present article derives the exact equations for the mean 
flow and the Reynolds stresses for a compressible turbulence keeping the neglected difference 
between the Reynolds and the Favre mean velocities. In so doing it is hoped to make two 
crucial points clear: 1) that the retention of the mass flux terms in the several places it 
appears in the mean momentum, mean energy and Reynolds stress equations is essential to 
the prediction of any nominally complex engineering flows and 2) that a careful distinction 
between Reynolds-averaged and Favre-averaged variables must be made to properly pose the 
near-wall problem and insure that the near-wall asymptotics are carried out consistently. 

In this exposition upper case letters will be used to denote mean quantities except in the 
case of the mean density, < p >, as p has no convenient upper case form. The expectations 
will be indicated using the angle brackets for time-means, eg. < >, and the curly brackets 

for the density-weighted or Favre-means, eg. {u.Wj}, where {u.u,} =< p*v±Vj >< p > -1 and 
the asterisk denotes the full field, p* =<p> +p r . The dependent variables are decomposed 


according to 

u* = Ui + Ui 

where < u, >= 0 


u* = Vi + Vi 

where {«;} = 0 


p* = < p> + p' 

where < p' >= 0 


p* = p + p 

where < p >= 0. 


As both the Reynolds-mean and the Favre-mean velocities appear in the evolution equations 
for a compressible turbulence it is necessary to carry both the Favre and the Reynolds 
decompositions of the velocity field. The mean of the fluctuating Favre velocity, < u, >, 
characterizes the difference between the Favre mean velocity and Reynolds mean velocities, 
Vi and Ui, as well as the difference between the instantaneous fluctuating velocities: 

Ui = v; + < Vi > 

U, = Vi — < Vi > . 

The Favre fluctuation mean, a first-order moment, is related to the turbulent mass flux, a 
second order moment, by 

< pvi >= — < p X Vi > 
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which follows from the definition of the Favre-average of the Favre fluctuation, {v,} =< 
p*Vi >< p > -1 = 0. Thus, apart from a scaling by the local mean density, the mean of the 
fluctuating Favre velocity and the mass flux are equivalent quantities. The primes on the 
fluctuating density have been dropped. 

Mathematically < u, > represents the difference between unweighted and density- weighted 
averages of the velocity field and is therefore a measure of the effects of compressibility 
through variations in density. It plays an important role in parameterizing the anisotropic 
effects of compressibility associated with the mean dilatation and the mean density gradi- 
ents. Experimentally it is an important and essential quantity that allows numerical results 
computed in Favre averaged variables to be related to experimental results computed in 
Reynolds variables. Additional results as well as an equation for the mass flux have been 
given in Ristorcelli (1993). The present investigation focuses on the role the Favre fluctuation 
mean plays in the the mean momentum, mean energy and the Reynolds stress equations and 
shows that the current practice of neglecting it is unnecessary, inconsistent and inadequate. 

The mean flow equations 

* . 

Substituting in the Favre decomposition into the Navier-Stokes equations and time av- 
eraging produces, without approximation, 

(< P > K),t + (</?> V p Vi), p = -P,i +£ tJ ,j ( U ) - (< p > {u,u P }), P 

where £„(t/) =< p > [Ui,j +U,,i — 2/3 Correlations between the fluctuating viscos- 
ity and the fluctuating velocity have been neglected. Note that the viscous terms are given 
in terms of the time-averaged mean velocity while the problem is solved in terms of the 
Favre- averaged mean velocity. They are related by (/,• = Vi + < u,- >. The usual assump- 
tion used to close the viscous term is that (/, ~ V{. Thus even the first-order equation is 
modeled reducing the accuracy of the method at the very lowest order in the very region of 
the most practical (aerodynamic) interest. Part of the appeal of k — e or Reynolds stress 
turbulence models in the incompressible turbulence is that the equations for the mean flow 
are carried exactly: this can also be done for the compressible turbulence. Substituting 
{/, = Vi + < Vi > produces the exact equation for the mean momentum equation: 

(< P > V<),« +(<P> V p Vi), p = -Pp + [Ey+ < <Tij >],; -(< p > {u,Up}),p 

where = <C p > k) 2/3 «md <C c xj > == <C p > [<C >,j T Vj > „■ 2/3 <C 

v g The mean momentum equation is now carried exactly. Note that the additional 

terms, < cr^ >, do not contribute to the mean flow equations if the turbulence is either 
homogeneous or isotropic. In the near-wall region where the viscous terms are important 
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and where there are gradients in the mean density, the < > are important and need 

to be carried. Data from Ma = 4.5 DNS computations of Dinavahi and Pruett (1993) in 
unidirectional developing wall bounded flow indicate that the approximation of f/ t ~ Vi in 
the wall bounded flow is inadequate. This is a nominally simple flow, in comparison to 
those of practical interest, in which there is a four-fold variation of the mean density over 
the boundary layer. In data taken from that simulation, shown in Figure 1, it was quite 
unexpectedly found that < v 2 > is larger than either U% and It is large enough to cause 
t/ 2 and V 2 to have different signs. This indicates that the mean fluid particle transfer is 
in a direction opposite to the net momentum transfer in flows with mean density gradients. 
Moreover, the contribution of the Favre fluctuation mean to the total viscous stress was found 
to be as large as one third the contribution of the Favre mean viscous stress, in the near 
wall portions of the turbulent boundary layer. In flows with separation and re-attachment 
the second derivative of the mean flow vanishes leaving, in the viscous terms, the second 
derivative of the Favre fluctuation mean. Clearly this will be an important term when there 
is a turbulent mass flux due to mean density gradients near the point of separation. 

The mean energy equation is subject to similar deficiencies. The exact equation for the 
mean total energy, internal plus kinetic, is 

(< p > E), t + (< p> V p E), v = - [PV p + P < v p > + < pv p >], p + [< VpkVk > - < p> {eu p }],p 

+ [SpfcVfc-j- < Gpk 14 Ep£ Vk >]>p Qv>v • 

It is not unusual in compressible turbulence models to see several, if not all, of the terms 
involving the mass flux, P < Vk >, < cr p k >14, or E p * < u* > dropped because the present 
models for these terms destabilize computations or because they are, in the spirit of t/, — K‘, 
assumed negligible. 

The Reynolds stress equations 

The second moment equations for a compressible flow, are written without approxima- 
tion, and after some manipulation, as 

D/Dt (< p > {uiUj}) = - < p > {viVp}Vj, p - < p > {vjV p }Vi, p + n.i +2/ 3 < P v k,k > tij 

-[< pVi > 8pj+ < pVj > 6 ip + <P> {viVjVp}- < Vj<T? p > - < Vi<Tj p >],p 

T ^ > [ P ii ~ b ^ikik "I" ^ &ik >jfc]"b ^ > [ f 5 ij T Ejkik "h ^ &jk >><: ] 

- < Uj, p (T > - < Ui, p <Tj p > 

where the mean momentum equations have been used and <r£ =< p, > [ui,j +Uj„- —2/3 u q , q 
The form of the equations above reflect the following manipulations: 1) The deviatoric part 
of the pressure-strain correlation is defined fltj p( v tij + u jji ) > —2/3 < pvk,k> & ij • an d 
2) the identity v, = u* + < u,- > has been used to rewrite the transport terms in Vi variables 
while keeping the dissipation terms in Ui variables. In the Reynolds stress equations the 
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terms arising from viscous surface terms appear naturally in u, variables while the problem 
is posed in the Favre v, variables. 

The dissipation type terms, which require closure, are kept in Ui variables. The fact 
that Vi = u, + < u, > is used to segregate the terms into ones that require modeling 
and ones that are carried as dependent variables in the closure scheme. The identities 
^ u^Up ^ ,ip 2 ^ Up ^ ,p 4" ^ u,',, - ^p?p ^ and Uj,p ^p^i where the 

vorticity is uij = tijkUj,k are used to write the trace of the dissipation terms as 

<C Ujip <7j*p !> =< C n ^ [<C Lv^ufk ~b 2(<C UjUp ,jp 2 <C u/-,^ Up ^>,p ) T 4/3 <! u p ^ p >]. 

Here <C DjL Oj ]> = < C 1ij,pUj,p^> Uj^pUp^j !> = <C Uj,pUj,p^> UqUp > ,qp ~i~2 <[ Uq,qUp > 

,p — < u g „Up,p>. In a homogeneous compressible turbulence < u^uj, >=< Uj, p Uj, p > — < 
Uj,j u p ,p > which reduces to the usual expression in an incompressible homogeneous turbu- 
lence. Defining the positive definite scalar dissipation quantities, the solenoidal dissipation 
< p > e a =< // >< WjWj > and the dilatational dissipation < p > Ed = 4/3 < // >< 
u p,pU g ,q > the trace can then be written as 

< Uj, p o/ p >= < p > (e 4 + £d ) + 

where e< = 2 < // > [< >,j —2 < u*,* u p >], p is a scalar transport term that comes from 

the dissipation type terms. The decrease of the Favre kinetic energy, for a homogeneous 
turbulence is then written as D/Di (< p > {vjVj}) = —2 < Uj, p o/ p >= —2 < p > 
(£s + £d) — 2e t where E t = 0 Note that e t is not defined as per unit mass quantity to emphasize 
that an equation need not be carried for it: the substitution u, = t>, — < u, > shows that 
it can be written, except for the correlation with the divergence, in terms of the mass flux 
and the Reynolds stresses for which equations are carried. In general Et is either positive 
or negative: in a homogeneous turbulence £t is zero while in the near-wall region it makes 
a nonnegligible contribution to the energy budget. The point is that in an inhomogeneous 
turbulence it is necessary to recognize the contributions to < Uj, p o/ p > that are functions of 
gradients of the Reynolds stresses rather than hiding these terms in the dissipation. If the 
trace is added and subtracted the dissipation-type terms can be rewritten as 

- < Uj, p o/ p > - < Ui, p >=,- < p > efj - 2/3 < p > (e, + Ed)Si } - 2/3 s t 6ij 

where the term < p > ef, has zero trace. The viscous terms are now manipulated into their 
final form. Using < 7 - — Oij - < cr,j >, which follows from u, = u, — < u, >, allows the 
viscous transport terms to be rewritten as < V{Oj p >=< ViOj p > — < Vi >< Oj p >. The 
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Reynolds stress equations can be written in the following form 
D/Dt (<p> {wit;,-}) = - < p > {viV p }Vj, p -<p> {vjV p }V itp + n*j +2/3 < pd > 6 {j 

' -[< pvi > 8 pj + < pvj > 8 ip + < P > {viVjV p}],p 
+ 2[< p >< vjVi >, p - < p >< Vid > 8 jp / 3- < p >< vjd > 8 iv j 3], p 
-2/3 e t 6ij - [< Vj >< <r ip > + < v,- >< (T jp >], p 

+ < Vj > [— P„ + S,p,p + < (Tip >,p] + < v i > [ — Pij + ^ ® jp >'P ] 

- < p > e d ij - 2/3 < p > (e, + e d )8ij 

where the definition of the stress tensor <Tij has been used to expand [< Vj(Ji p > + < Vi<Tj p >] 
in the viscous transport terms and d = v p , p . Except for neglecting correlations between the 
fluctuating Favre velocity and its vorticity the final form of the Reynolds stress equations, 
given above, is exact. 

Conclusions 

The mass fluxes or Favre fluctuation means appears in several places in the compressible 
turbulence equations. They contribute to 1) the viscous terms in the mean momentum 
equation, and in the mean energy equation they contribute to 2) the viscous terms, 3) the 
pressure work terms, 4) the viscous work terms, and the 5) pressure flux which is coupled to 
the mass flux through the equation of state. Modeling Ui — Vi is equivalent to neglecting the 
mass flux in five different locations in the mean equations. In the Reynolds stress equations 
the mass flux contributes to 6) the viscous diffusion of the Reynolds stresses, which only 
recognized when the viscous terms are properly distinguished into their Reynolds and Favre 
variable components. The mass flux determine the importance of two production mechanisms 
7) one due the acceleration of the mean flow and 8) the other due to viscous effects associated 
with the Favre fluctuation mean and also it contributes to 9) the pressure flux. A general 
algebraic model for the < Vi >, derived from first principles, and suitable for complex flows 
of engineering interest has been derived in Ristorcelli (1993) and tested in Ristorcelli and 
Dinavahi (1993). 

The viscous terms in the Favre-averaged Reynolds stress equations have been system- 
atically derived, identified and segregated. The source of terms associated with anisotropy 
and inhomogeneity of the dissipation, essential to consistent near-wall modeling, have been 
shown. The second essential point is that the dissipation is a function of u 4 fluctuations 
about the Reynolds mean, while the Reynolds stress equation is an evolution equation for 
statistics of fluctuations around the Favre mean. The near-wall Taylor series expansions 
of these two quantities, in a general compressible flow, are different. Near-wall asymptotics 
must recognize these facts when relating Favre type variables {u.Uj}, in which the prob- 
lem is posed, to dissipation type quantities which are carried in Reynolds variables. This 
point is crucial for any consistent general near-wall model development in flows only nomi- 
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nally more complex than the usual uni-direction adiabatic/isothermal flat plate flows: it is 
necessary to start with the exact equations distinguishing between Favre and Reynolds de- 
composed variables and recognizing that their respective Taylor series expansions are related 
by Vj = Uj-f < vj >. 

It is typical of current compressible k—e and Reynolds stress models to neglect some or all 
of these contributions to the mean momentum, mean energy and Reynolds stress equations. 
The retention of the mass flux terms will be important in complex compressible turbulent 
flows: these include flows in which there are mean density gradients due to large Mach 
number or combustion, separation or reattachment (inflection points), cold wall boundary 
conditions, mean dilatation, shocks, adverse pressure gradients, or strong streamwise accel- 
erations. 
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